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Abstract 

We study integrals of motion and factorizable S-matrices in two-dimensional integrable 
field theory with boundary. We propose the "boundary cross-unitarity equation" which is 
the boundary analog of the cross-symmetry condition of the "bulk" S-matrix. We derive 
the boundary S-matrices for the Ising field theory with boundary magnetic field and for 
the boundary sine-Gordon model. 



1. INTRODUCTION 

In this paper we study two-dimensional integrable field theory with boundary. Exact 
solution to such a field theory could provide better understanding of boundary-related 
phenomena in statistical systems near criticality[l]. Quantum field theory with boundary 
can be applied to study quantum systems with dissipative forces [2]. From a more general 
point of view, studying integrable models could throw some light on the structure of the 
"space of boundary interactions" , the object of primary significance in open string field 
theory [3]. 

An integrable field theory posesses an infinite set of mutually commutative integrals 
of motion 1 . In the "bulk theory" (i.e. without a boundary) these integrals of motion 
follow from the continuity equations for an infinite set of local currents. These currents 
can be shown to exist in many 2D quantum field theories, both the ones defined in terms 
of an action functional (like sine-Gordon or nonlinear sigma-models[4]) and those which 
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In practice, it is usually enough to find the first few integrals of motion of spin s > 1 to 

argue about integrability (this argument works in all known cases). 
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are defined as "perturbed conformal field theories" [5] . In presence of a boundary, the 
existense of these currents is not sufficient to ensure integrability. Integrals of motion 
appear only if particular "integrable" boundary conditions are chosen. In general, the 
boundary condition can be specified either through the "boundary action functional" or 
as the "perturbed conformal boundary condition" 2 . In Sect. 2 we show how the integrals 
of motion of the "bulk" theory get modified (in most cases destroyed) in presence of the 
boundary and how one can find "integrable" boundary conditions. 

An important characteristic of an integrable field theory is its factorizable S-matrix. 
In the "bulk theory", the factorizable S-matrix is completely determined in terms of the 
two-particle scattering amplitudes, the latter being required to satisfy the Yang-Baxter 
equation (also known as the "factorizability condition"), in addition to the standard equa- 
tions of unitarity and crossing symmetry [4]. These equations have much restrictive power, 
determining the S-matrix up to the so-called "CDD ambiguity" . At present many examples 
of factorizable scattering theory are known (see e.g. [4,9-12]), most of which are obtained by 
explicitely solving the above equations (eliminating the "CDD ambiguity" usually requires 
a lot of guesswork). 

It is known since long[13] that the concept of factorizable scattering can be generalized 
in a rather straightforward way to the case where a reflecting boundary is present. The 
S-matrix is expressed in terms of the "bulk" two-particle S-matrix and specific "boundary 
reflection" amplitudes, the latter, again, being required to satisfy an appropriate gen- 
eralization of the Yang-Baxter equation (which we call here the "Boundary Yang-Baxter 
equation") 3 . Generalization of the unitarity condition is also fairly straightforward. What 
was not known was the appropriate analog of the crossing-symmetry equation. In Section 3 
we fill this gap by deriving what we call the "boundary cross-unitarity equation" . Together 
with this, the above equations have exactly the same restrictive power as the corresponding 
"bulk" system, i.e. they allow one to pin down the factorizable boundary S-matrix up to 
the "CDD factors" . 

We also study integrable boundary conditions in two particular models (off-critical 
Ising field theory and sine-Gordon theory) and find the associated boundary S-matrices. 
This is done in Sections 4 and 5. 



2. INTEGRALS OF MOTION 

Consider a 2D Euclidean field theory, in flat space with coordinates (x 1 ,^ 2 ) = (x,y). 
There are basically two ways to define a 2D field theory. In the Lagrangian approach one 
specifies the action 



2 Conformal field theories with boundary are studied in [6,7,8] where, in particular, it is 
shown how to classify conformally-invariant boundary conditions and boundary operators. 

3 This equation found important applications in quantum inverse scattering method [14], 
generalized to the systems with boundary [15, 16]. 
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where <p(x, y) is some set of " fundamental fields" and the action density a(cp, d^ip) is a local 
function of these fields and derivatives d^ip = dip/dx^ with fj, = 1,2. Another approach is 
to consider the "perturbed conformal field theory"; in this case one writes the "symbolic 
action" 

/oo 
§(x,y)dxdy (2.2) 
-oo 

where Aqft is the "action of conformal field theory (CFT)" and <&(x, y) is a specific 
relevant field of this CFT. In both approaches one can define a symmetric stress tensor 
T M!/ = T Ufl which satisfies the continuity equations 



d- z T = d z Q ; d z T = d- z Q (2.3) 

Here we use complex coordinates z = x + iy,z = x — iy and denote the appropriate compo- 
nents T = T zz , T = T zz , © = T zz of the stress tensor. To achieve Hamiltonian formulation 
one chooses an arbitrary direction, say the y-direction, to be the "euclidean time" , and as- 
sociates a Hilbert space H with any "equal time section" y = const., x G (— oo, oo). States 
are vectors in H and their "time evolution" is described by the Hamiltonian operator 

/oo poo 
dxT yy = / dx[T + f + 26] (2.4) 
-oo J — CO 

Let us assume that the field theory (2.1) or (2.2) is integrable. In particular, the 
equations (2.3) appear to be the first representatives of an infinite sequence 



d z T s+1 = d z e s -t d- z f s+1 = a z e s _i (2.5) 

where T s+ i,© s _i (T s+ i, s _i) are local fields of spins s + l,s — 1 respectively and the 
integrals of motion (IM) 

/oo />oo 
(T s+1 + e s _i)rfx; P s = (f s+1 + 6 s _i) dx (2.6) 

-oo J — oo 

constitute an infinite set of mutually commutative operators in H. The spin s of IM (2.6) 
takes integer values s±, S2, ■■■ in the infinite set {s} which is an important characteristic of 
an integrable field theory [5]. In any case, s\ = 1; for this value of s (2.5) coincides with 
(2.3) and 



H=(P l + P 1 ) (2.7) 

Now, let us consider this field theory in the semi-infinite plane, x G (— oo,0], y G 
(—00,00), the y-axis being the boundary. Again, the boundary conditions are specified 
in different ways in the two approaches, (2.1) and (2.2). In the lagrangian approach 
one chooses the "boundary action density" 6(^e(y), ^<^s(y)), as a local function of the 
"boundary field" <pb, <fB(y) = ^p(x, y)\ x= o, and writes the full action in the form 
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/OO rO poo ^ 

dy dxa(ip,df,ip) + dy b(ip B , -t^b) (2.8) 
-oo J —oo J — oo dy 

To write down the analog of (2.2) in presence of the boundary, one starts with conformal 
field theory on the same semi-infinite plane with certain conformal boundary conditions 
(CBC) at the boundary x = 0, and defines "CBC perturbed by relevant boundary operator 
$s(y)". In general, this perturbation of boundary condition goes along with the pertur- 
bation (2.2) of the bulk theory. This strategy is summarised by the symbolic "action" 

/oo i>0 poo 

dy dx$(x,y)+ dy$ B (y) (2.9) 
-oo J — oo J — oo 

As is argued by Cardy[6], in CFT the equation T xy \ x= o = is satisfied with any choice 
of CBC. In the perturbed theory (2.9) this condition is changed to 

T xy \ x = = (-i)(T-f)\ x=0 = |-%) (2.10) 

where 9(y) is some local boundary field. As the theory (2.8) or (2.9) is still symmetric 
with respect to translations along the y-axis, the equation (2.10) can be easily derived as 
a consequence of this symmetry. In the "perturbed CFT" approach, it is relatively easy 
to relate the field 9{y) to the "boundary perturbation" $ B (y)(see below). 
The continuity equations (2.3) guarantee that the contour integrals 

P 1 {C)= [{Tdz + Qdz) ; P^C) = I \fdz + Qdz) (2.11) 
Jc Jc 

do not change under deformations of the integration contours C. Consider the contour 
C = Ci + C12 + C 2 shown in Fig.l. Obviously, Pi(C) = A(C) = 0. If we take the 
combination 



= Px(C) + Pi(C) = Pi(d) + P!(Ci) + Pi(C 2 ) + Pi(C 2 ) + P!(C 12 ) + P!(C 12 ) (2.12) 
the integration over C 12 part of this contour is easily done in view of (2.10) 

Pi(C i2 ) + Pi(C 12 ) = 9( yi ) - 9(y 2 ) (2.13) 

and hence the integral 

H B {y)= f (T + T + 2Q)dx + 9{y) (2.14) 

J — oo 

is in fact y-independent, i.e. it is an integral of motion in (2.8) or (2.9). 

Even if the bulk theory (2.1) or (2.2) is integrable, in general the boundary conditions 
in the semi-infinite system (2.8) or (2.9) will spoil integrability. Suppose, however, that 
we can choose particular boundary conditions, such that the equation 
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[t 8+1 + e s _! - f s+1 - e s _i]U =0 = -j-Osiy) (2.i5) 

is satisfied for some s G {s}; here again 9 s (y) is some local boundary field. Then repeating 
the above argument, one finds that the quantity 

H%> = [T 8+1 (x,y) + e s - 1 (x,y) + f s+1 (x,y) + d s - 1 (x,y)}dx + 9 s (y) (2.16) 

J — oo 

does not depend on y, i.e. it appears as a non-trivial integral of motion. We will call 
the boundary conditions "integrable" (and refer to the theory (2.8) ((2.9)) as "integrable 
boundary field theory") if the equation (2.15) holds for any s out of an infinite set {s}b, 
subset in {s}. 

There are two alternative natural ways to introduce the Hamiltonian picture in the 
theory (2.8)((2.9)). First, one can take again the direction along the boundary (y-direction) 
to be the "time". In this case the boundary appears as the "boundary in space", and the 
Hilbert space of states H.b is associated with the semi-infinite line y = const, x G (— oo, 0]. 
Then the quantities (2.14), (2.16) appear as operators acting in H,b, and Hg(= 
is naturally identified with the Hamiltonian. The correlation functions of any local fields 
Oi(x, y) in presence of the boundary can be computed in this picture as the matrix elements 



/n f \ n ( \\ B \ 7 y0 1 (x 1 ,y 1 )...0 N (x N ,y N ) \0) B 
(0 1 (x 1 ,y 1 )...0 N (x N ,y N )) = y - — — (2.17) 

B\U I V/B 

where | 0)sG TCb is the ground state of Hb, and Oi(xi,yi) in the r.h.s. are understood as 
the corresponding Heisenburg operators 

Oi(x, y) = e- yHB O l {x, 0)e yHB (2.18) 

and T y means the "y-ordering" . In an integrable theory the operators H^; s G 
constitute a commutative set of IMs. 

Alternatively, one could take x to be the "euclidean time". In this case the "equal 
time section" is the infinite line x = const, y G (—00,00). Hence the associated space of 
states is the same Ti as in the bulk theory (2.1)((2.2)), and the Hamiltonian operator is 
given by the same eq.(2.4) (with x and y interchanged). The boundary at x = appears as 
the "time boundary" , or "initial condition" at x = which is described by the particular 
"boundary state" 4 ] B) G H. It is the state | B) that concentrates all information about 
the boundary condition in this picture. The correlation functions (2.17) are expressed as 

in i \ n ( \\ (0 I %(0 1 (x 1 ,y 1 )...0 N (x N ,y N ) | B) 
{0 1 (x 1 ,yi)...0 N {x N ,y N )) = (2.19) 



4 The notion of boundary state is discussed in the context of CFT in [7] . 
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where now | 0) G H is the ground state of H and Oi(x, y) in the r.h.s. are the Heisenberg 
field operators 

O i (x,y) = e- xH O i (0,y)e xH (2.20) 

corresponding to this picture; T x means "x-ordering" . In an integrable theory, the same 
equations (2.6) (again with x and y interchanged) define an infinite set of mutually com- 
mutative operators P s , P s ; s G {s}, acting in H. As a direct consequence of (2.15) one finds 
that the boundary state | B) satisfies the equations 

(P s -P s )\B) = ; se{s} B . (2.21) 

To expose an example of an integrable boundary field theory, let us cosider the per- 
turbed CFT (2.9) with Ac ft taken to be any c < 1 minimal model, and the degenerate 
spinless field $(1,3) taken as the bulk perturbation, 

*(x,y) = A* (1)3 )(x,y). (2.22) 

where A is a constant of dimension [length] ; A = A( 13 ). The local integrals of motion 
in the corresponding bulk theory (2.2) are discussed in [5]. The fields T s+ i are composite 
fields built up from T = T zz and its derivatives, for example 

T 2 =T ; T 4 =:T 2 : ; T 6 =: T 3 : : (d z T) 2 : ; ... (2.23) 

where : : denotes appropriately regularized products (T s+ i are built from T in similar 
ways). The characteristic feature of the fields T s+ i is that, in the conformal limit A = 
their OPE's with $(1,3) have the form 

s \l/ (fc) (w w) 

T s+1 {z)$ (li3) {w,w) = J2 a 7 _ \k + regular terms (2.24) 

fc=i ^ W > 

(^s-fe+i are particular conformal descendants of $(1,3)) with 

^\w,w) = d w Q s - 1 {w,w), (2.25) 
where Q s -i are local fields. The "null- vector" equation 

(£-3 - ^-iL-* + (A + 1 ) (A + 2) ^.)*W < 2 ' 26 ) 

satisfied by the degenerate field $(1,3), is used to prove (2.25). The equation (2.25) is 
sufficient to show that the field T s+ i satisfy (2.5) in the perturbed theory, up to first order 
in A. Then dimensional analysis shows that (2.5) is exact. There are infinitely many fields 
T s+ i satisfying (2.24)-(2.25), one for each odd s. So, the set {s} in this theory contains all 
positive odd integers. 

Conformal boundary conditions in c = 1 — p ^ +1 ^ minimal CFT are classified in [7]. 
There are finitely many possible conformal boundary conditions, each corresponding to 
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a particular cell (r, s) of the Kac table. In each case possible boundary operators are 
degenerate primary boundary fields ^(n.m) (plus their Virasoro descendants), such that 

the fusion rule coefficients -^(„^)( rs ) are non-zero (the other degenerate boundary fields, 

with ^nm)(r s) = ® correspond to "juxtapositions" of different conformal boundaries, see 
[7] for details). The field $(1,3) satisfies this condition for any (r, s). We take this field to 
be the boundary perturbation in (2.9), i.e. 

= A B V(i,3)Q/)- (2-27) 

Its conformal dimension A = A( 13 ) < 1, so that it is a relevant perturbation. We want to 
show that under this choice the theory (2.9) is integrable. 

To warm up, let us consider the components T and T of the stress tensor itself. In 
the conformal limit A = 0, Xb = 0, these components satisfy the boundary condition 

[T(y + ix) - f(y - ix)]\ x=0 = 0, (2.28) 

i.e. the field T(z) is just the analytic continuation of T(z) to the lower half-plane Imz < 0. 
Let us "turn on" the boundary perturbation, A^^O, still keeping A = 0, and consider the 
correlation function 

([T(y + ix)-T{y-ix)]X) XB = 

Z^([T(y + ix) - T(y - lX )]Xe- XB *™«W) CFT , (2.29) 
where X is any product of fields located away from the boundary x = and = 

(e Xb ^ (1 ' 3) ^ dy ) CFT . In the limit x — > the contribution to (2.29) is controlled by 
OPE 

(T(y + ix) - T(y - ix))\ B ip(i,3)(y') = 

Xb r A (i,3) A (i,3) + 1 g 1 9 ^ ^ , ^ 

B {y — y' + ix) 2 (y — y' — ix) (y — y' + ix) dy' (y — y' — ix) dy' 

- X B {A (h3) 6'(y - y') + 5(y - y ')-^}^ 3) ( y '). (2.30) 
This shows that for Xb ^ the fields T(x, y) and T(x, y) satisfy (2.10) with 

%) = (l-A)A B V(i,3)(</). (2-31) 

Now, dimensional analysis, exactly parallel to that carried out in [5] shows that the equa- 
tions (2. 10), (2. 30) remain valid if we turn on the bulk perturbation, i.e. at A 7^ 0. 

The above arguments can be repeated for the higher currents T s +i, s = 3, 5, .... To see 
this, note that the null vector equation (2.26), crucial for the validity of (2.25), is satisfied 
by ^(1,3) as wen - Therefore, in the conformal limit A = Xb = the fields T s+ i(z) satisfy 
the OPEs 
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s 1 

T a+ i(z)V>(i, 3 )(y) = 5^ , _ N fc xi+i-fc(y) + regular terms, 
fc=i ^ yJ 

xi 1) (y) = ^ s (y), (2-32) 

similar to (2.24), (2.25) where are descendants of ^(1,3) • Whence 

\imit x ^ (T s+1 (y + ix) - T s+1 (y - ix)) \b^{i,3){v') 
rJ . s - 1 

= \ B (5(y - v')^Mv') + H {k - 1)ldyt k-A v - y')x k s+1 - k (y% (2.33) 

and we conclude that for A = and in the first order in A# the equation holds 

[T a+1 (x,y)-T a+1 (x,y)] x=0 = ^9 8 (2.34) 

with 

s (y) = X B [ qs (y) + E (fc -i)i^-i xi + i-*(y)]- (2-35) 

The higher powers of As can contribute to (2.34) through the "resonance terms" similar 
to those discussed in [5]. It is plausible, however, that these do not spoil the general form 
of (2.34) but simply modify (2.35) by higher order terms in As. It is also plausible (and 
can be supported to some extent by dimensional analysis of [5]) that "turning on" the bulk 
perturbation, A ^ 0, converts (2.34) to(2.15). 

We realise that the above arguments do not constitute a rigorous proof. First, we 
did not solve the problem of the "resonance terms" (this problem remains open in the 
bulk theory, too). More importantly, we did not analyse possible effects of mixing between 
the boundary and the bulk perturbations. We have shown,however,that the field theory 
(2.9) with (2.22) and (2.26) satisfies some very non-trivial necessary (but not sufficient) 
conditions of integrability and we conjecture that this boundary field theory with boundary 
is integrable. 



3. BOUNDARY S-MATRIX 

If the field theory (2.1) ((2.2)) is massive the space H is the Fock space of multiparticle 
states. After rotation y = it to 1 + 1 Minkowski space-time these states are interpreted 
as the asymptotic ("in-" or "out-") scattering states. For an integrable field theory the 
scattering is purely elastic and the corresponding S-matrix is factorizable. The factorizable 
scattering theory in infinite space is discussed in many papers and reviews (see e.g. [4,5, 9- 
12]); below, we describe just some basics of the theory. The boundary theory (2.8) ((2.9)) 
in Minkowski space is also interpreted as a scattering theory. For the integrable boundary 
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field theory this scattering theory is again purely elastic and the corresponding S-matrix 
is the "factorizable boundary S-matrix". The "factorizable boundary scattering theory" 
is developed in close parallel with the "bulk" theory (see [13]). However there are still 
some gaps in this parallel (most importantly, the boundary analog of crossing-symmetry 
condition of the "bulk" S-matrix is not absolutely straightforward). It is the aim of this 
Section to fill these gaps. 

We start with a brief description of the basics of the factorizable scattering theory in 
the infinite space ("bulk theory"). Assume that the theory contains n sorts of particles 
A a ; a = 1,2, ...,n with the masses m a . As usual, we describe the kinematic states of the 
particles in terms of their rapidities 0, 

Po+Pi=me 9 ; p - p 1 = me~ e , (3.1) 

where p^ are the components of the two-momentum and m is the particle mass. The 
asymtotic particle states are generated by the "particle creation operators" A a (0) 

I A ai (e 1 )A a2 {e 2 )...A aN (e N )) = A ai (e 1 )A a2 (9 2 )...A aN (e N ) \ o>. (3.2) 

The state (3.2) is interpreted as an "in-state" if the rapidities 9i are ordered as 9\ > 2 > 
... > 9n', if instead 9\ < 2 < ••• < On (3.2) is understood as an "out-state" of scattering. 
The "creation operators" A a {9) satisfy the commutation relations 

AaMAaM = S^fr - 9 2 )A b2 (9 2 )A bl (9 1 ) (3.3) 

which are used to relate the "in-" and the "out-" bases and hence completely describe the 
S-matrix. The coefficient functions (9) are interpreted as the two-particle scattering 
amplitudes describing the processes A ai A a2 — > A 0l A 02 (see Fig. 2). The asymptotic states 
(3.2) diagonalise the local IM (2.6), the eigenvalues being determined by the relations 

[P 8 , A a {0)\ = la ^e se A a (0); [Ps, MO)] = l^e^ A a (9); (3.4) 

P s I 0) = 0; P s I 0) = 0, (3.5) 

where 7a are constants (7! = n^a)- These IM must commute with the S-matrix; it 
follows in particular that the amplitude S%\a (0) is zero unless m ai = rrib 1 and m a2 = 
(other concequences are discussed in [5]). 

Charge conjugation C acts as an involution of the set of particles {A a }, i.e. C : A a <-> 
Aa, where A a E{A a }, so that each particle in {A a } is either neutral CA a = A a or belongs 
to the particle-antiparticle pair (A a , A a ). In this Section we assume for simplicity that the 
theory under consideration respects C, P and T symmetries, i.e. 

S b a \ b a \ (0) = S\fi 2 (0) = S b a f a \ (0) = Sgf (0). (3.6) 

The two-particle S-matrix S^\^ (6) is the basic object of the theory. It must satisfy 
several general requirements 

1. Yang-Baxter (or "factorization") equation 
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si\ c a\{e)s b c iz{e + e')s b c f c i{9') = szz(9')s c a ii(9 + 9')s b c \ b i(9)- (3.7) 

here and below summation over repeated indices is assumed. This equation is illustrated 
in Fig. 3. Formally, this equation appears as the associativity condition for the algebra 
(3.3). 

2. Unitarity condition 

^ C a 2 2 W^ 1 1 c 2 2 (-^)=^i^ 2 - (3-8) 

Graphic representation of this equation is shown in Fig. 4. It can also be obtained as the 
consistancy condition for the algebra (3.3) (one applies (3.3) twice). 

3. Analyticity and Crossing symmetry. The amplitudes S blb ? (6) are meromorphic 
functions of 9, real at Re9 = 0. The domain < Im9 < re is called the "physical strip". 
The physical scattering amplitudes of the "direct channel" A ai A a2 — > A^A^ are given by 
the values of the functions S b a lb ^ 2 (9) at Im9 = 0, Re9 > 0. The values of these functions at 
Im9 = 0, Re9 < describe the amplitudes of the "cross-channel" A a2 Aj Ji — > A h2 Aa x . The 
functions S blb a 2 2 (9) satisfy the crossing symmetry relation 

Slim = Slf h [{™ - 9) ( 3 - 9 ) 
(see Fig. 5). Combining this and the Eq.(3.8) one can derive the "cross- unitarity equation" 

S C a ic 2 M ~ ^atS™ + *) = (3-10) 

4. Bootstrap condition. The only singularities of S blb ^ 2 (9) admitted in the physical 
strip are poles located at Re9 = 0. The simple poles are interpreted as bound states, 
either of the direct or of the cross channel. As the bound states are stable particles they 
must be in the set {A a }. Let iu c aia2 be the position of the pole of S^ b ^ 2 {9) associated with 
the "bound state" A c of the direct channel. Then u c aiCL2 must satisfy the relation 

TO a x + ™l 2 -™ 2 c = -2m ai m a2 cos< ia2 , (3.11) 

i.e. the quantity u c aiCL2 = n — u c aia2 can be interpreted as the internal angle of the euclidean 
triangle with the sides m ai , m a2 , m c . The pole term 

fc fb 1 b 2 

s b ai\{e)^i J ;^ Jc c (3.i2) 

corresponds to the diagram in Fig. 6, where the vertices represent the "three-particle cou- 
plings" / (Fig. 7). In this situation the two-particle S-matrix satisfies the "bootstrap equa- 
tion" 

f C a ia2 S bb m = flc 2 S c a i\{9 + iull,)SZ%{0 - i«£ e ) (3.13) 

which is illustrated by Fig. 8. 

More details about the factorizable scattering theory and many examples can be found 
in the original papers and reviews (see e.g. [4,5,9-12]). Most of the examples are constructed 



10 



directly, by solving the Eq. (3.7-3.9) above. Following this approach one can pin down the 
S-matrix S h a ^{6) up to the so-called "CDD ambiguity" 

(3-14) 

where the "CDD factor" <E>(0) is an arbitrary function satisfying the equations 

$(0) = $(in-e); $(#)$(-#) = 1. (3.15) 

The bootstrap equation may impose further restrictions on this function. 

Let us turn now to the semi-infinite system (2.8) ((2.9)). Here again the states in Hb 
can be classified as asymptotic scattering states. The scattering occurs in the semi-infinite 
1 + 1 Minkowski space-time (x, t),t= —iy, x < 0. The initial state 

| A ai (6 1 )A a2 (6 2 )...A aN (6 N )) B , in (3.16) 

(the subscript B indicates that | ...)# G Hb) of the scattering consists of some number 
(N) of "incoming" particles moving towards the boundary at x = 0, i.e. all the rapidities 
6>i, 6>2, 6>tv are positive. In the infinite future, t — > oo, this state becomes a superposition 
of the "out-states" 

I A bl (e 1 ')A b2 (e 2 ')...A bM (e M '))B,out (3.17) 

each containing some number of "outgoing" particles moving away from the boundary 
with negative rapidities 6\ , 62 ', 0m' '• In integrable boundary field theory this process is 
constrained by the IM (2.16). Like in the "bulk" theory the operators H s are diagonal in 
the basis of asymptotic states and 

H s I A ai (9 l )A a2 (6 2 )...A aN (6 N )) 

B ,in(out) 

N 

(X> 7 <>> cosh(^) + h^) 1 A ai {e 1 )A a2 {e 2 )...A aN {e N )) B ^ n{out) (3.18) 
i=i 

where are some constants. The constraints 

N M 

^ 7 (/cosh(^) = ^ 7 Wcosh«) (3.19) 
i=i i=i 

which follow from (3.18) show that M = N and the set of rapidities {(V, 6*2', 6* n'} can 
differ only by permutation from {— 61, —6*2, — 0/v}, i.e. the boundary scattering theory 
is purely elastic. It is possible to argue that the S-matrix in this case has a factorizable 
structure. 

The factorizable boundary scattering theory can be described in complete analogy 
with the "bulk" scattering theory. Again, the asymptotic states (3. 16), (3. 17) are generated 
by the "creation operators" A a {9) satisfying the same commutation relations (3.3). If 
B\ > 6 2 > ••• > On > the in-state (3.16) can be written as 
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A ai {e 1 )A a2 {e 2 )...A aN {e N ) \ o) B , (3.20) 

where | 0}s is the ground state of H B . One can think of the boundary as an infinitely 
heavy impenetrable particle B sitting at x = and formally write the state | 0)s as 

\0)b = B\0) (3.21) 

in terms of "operator" B which we call the "boundary creating operator" (formally B : 
7i — > Hb', it is an interesting question wether (3.21) makes any more than just formal 
sense). The "operator" B satisfies the relations 

A a (e)B = R b a (9)A b (-9)B, (3.22) 

the coefficient functions R a {9) being interpreted as the amplitudes of one-particle reflection 
off the boundary, as shown in Fig. 9. The Eq.(3.18) is reproduced if we assume (3.4), (3.5) 
and [H g ,B] = h {s) B. It follows from (3.17) that R h a {6) vanishes if m a ^ m b . By purely 
algebraic manipulations, with the use of relations (3.3) and (3.22), one can expand any 
in-state (3.20) in terms of the out-states 

A bl (-e 1 )A b2 (-e 2 )...A bN (-e N ) I o) B , (3.23) 

(9 1 > 9 2 > ... > On > 0) thus expressing the iV-particle S-matrix 

A ai (9 1 )A a2 (9 2 )...A aN (9 N ) \0) B = 

I 0)b (3-24) 

in terms of the "fundamental amplitudes" S^ 1 i b a 2 2 (9) and R b a (9) which are basic objects of 
the factorizable boundary scattering theory. The amplitudes R b a {9) have to satisfy several 
general requirements analogous to the requirements 1 — 4 of the "bulk" theory above. 
V . Boundary Yang-Baxter equation 

Kl^S:^ + 9 2 )Ril(9 1 )S b d f d \(9 l - 9 2 ) = 

S c a \ c a \ (9! - 9 2 )R d c l {9 1 )S d c l b l (9 1 + 9 2 )R% (9 2 ) (3.25) 

(represented graphically in Fig. 10) can be obtained as the associativiy condition of the alge- 
bra (3.22), (3.3). These equations have been introduced first in [13] and studied in relation 
with the quantum inverse scattering method for the integrable systems with boundary in 
many subsequent papers (see e.g. [15-17]). Note that (3.25) is the direct analog of (3.7). 
2'. Boundary Unitarity condition 

K{0)R b c {-9) = 5 b (3.26) 

is also an absolutely straightforward generalization of (3.8) (see Fig. 11). One obtains (3.26) 
applying (3.22) twice. 
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3'. The boundary analog of crossing-symmetry condition (3.9) is far less straightfor- 
ward. Note that without any additional conditions the equations (3.25) and (3.26) are not 
restrictive enough. The ambiguity in the solution is 

R h M -> R b a (9)$ B (9) (3.27) 
with arbitrary function $s which satisfy the equation 

$b(9)$b(-9) = 1 (3.28) 

only, which does not even imply (contrary to (3.15)) that $s is analytic (meromorphic) 
in the full complex plane of 9. 

To reveal the analog of the "cross channel" of the scattering process (3.24) one has to 
use the alternative Hamiltonian picture for the boundary field theory mentioned in Sect. 2. 
Consider 1 + 1 Minkowski space-time (r, y) with r = ix(r > 0) interpreted as time. The 
"equal time" section now is the infinite line — oo < y < oo and the space of states H is the 
same as in the "bulk" theory. The boundary condition at x = appears in this picture as 
the initial condition at r = 0; it is described by the "boundary state" | B) as explained 
in Sect. 2. As | B) e 7i this state is a superposition of the asymptotic states (3.2) of the 
"bulk" theory. In integrable theory with "integrable boundary" the states (3.2) admitted 
to contribute to | B) are restricted to satisfy (2.21). Noting that the eigenvalue of the 
operator P s — P s on (3.2) is 

N 

£27<M sinh(^) (3.29) 
i=i 

we conclude that the particles A a can enter the state | B) only in pairs A a (9)Ab(— 9) of 
equal-mass particles of the opposite rapidities 5 . Thus we can write 



\B)=MY, / d9 l d9 2 ...d9 N K^ aN ~ 1 - aiMb2 - bN (9 1 ,9 2 ,...,9 N ) 

A aN (-9 N )A aN _ 1 (-9 N - 1 )...A ai (-9 1 )A bl (9 1 )A b2 (9 2 )...A bN (9 N ) | 0) (3.30) 

where we have chosen the expansion in terms of the out-states; K 2 n are certain coefficient 
functions which can be related to the amplitudes R in (3.24). The overall factor Af in (3.30) 
is chosen in such a way that Kq = 1; it is possible to argue that in the massive theory 
with a non-degenerate ground state | 0)s, which we consider here, one can choose Af = 1 
by adding an appropriate constant term to the boundary action density b in (2.8) (or to 
the "perturbing field" in (2.9)); in what follows we assume this choice. By applying 
the standard reduction technique to the equality (2.17), (2.19) one can show that (under 
appropriate normalization of the operators A a (9)) the following equations hold 

5 The possibility of zero rapidity particles A a (0) entering the boundary state, which is 
evidently consistent with (2.21), is discussed below. 
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K a N a N . 1 ...a 1 , bl b 2 ...b N{du ^ ^ ^ = Rfrfc^gL _ Q u ™ _ l JL _ (3 . 31) 

up to a constant phase factor which depends on the normalizations of the operators A a (9); 
in what follows we assume that they are normalized in such a way that (3.31) holds as it 
stands. 

Let us concentrate attention on the amplitude K ab (9) = K2 ,b (9) describing the con- 
tribution 

poo 

\B) = (1 + K ab (9)A a (-9)A b (9) + ...) | 0) (3.32) 
Jo 

of the two-particle out-state A a (—9)A b (9) | 0) to the boundary state | B). It satisfies 

K ab (9)=R b d C{-0) (3.33) 

i.e. the "elementary reflection amplitude" R\ (9) can be obtained by analytic continuation of 
the amplitude K ab (9) to the domain ImO = if, Re9 < 0. The values of K ab {9) at negative 
real 9 (corresponding to the "lower edge" of the cut in the energy plane) are interpreted 
as the coefficients of expansion of | B) in terms of the in— states A a (9)A b (—9) | 0), 9 > 

/•OO 

\B) = (1+ / K ab {-9)A a {9)M~9) + •••) I 0) (3.34) 
Jo 

As the in— and the out— states are related through the S- matrix, the amplitude K ab (9) 
has to satisfy the following "boundary cross-unitarity condition" 

K ab {9) = S% (29)K b ' a ' (-9) (3.35) 

which we consider to be the boundary analog of (3.10). It is illustrated by the diagram 
in Fig. 12. With this equation added, the ambiguity in the solution of (3.25), (3.26) and 
(3.35) reduces to (3.27) with $ B satisfying (3.28) and 

® B ($) = $ B (in-9) (3.36) 

which is exactly the same as the "CDD umbiguity" (3.14), (3.15) in the "bulk" theory. 
Let us note that (3.35) allows one to write (3.32) as 

1 f°° 

\B) = (1 + - K ab (9)A a (-9)A b (9) + ...) | 0) (3.37) 

^ J-oo 

Using the Equation (3.31) one can express all the amplitudes K 2 n in (3.30) in terms of 
the two-particle boundary amplitudes K ab {9) and the elements of the two-particle S-matrix 
S^f(9). The result is in exact agreement with the following simple expression 

/OO 
d9K{9)) | 0) (3.38) 
-oo 
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where 



K{9) = -K ab (9)A a (-9)A b (9) (3.39) 

Note that although the "creation operators" A a (9) do not comute, the bilinear expressions 
(3.39) satisfy the commutativity conditions 

[K(9), K(9')] = (3.39) 

as a direct consequence of the "boundary Yang-Baxter equation" (3.25) and (3.33); there- 
fore there is no ordering problem in (3.38). The commutativity (3.39) is crucial for the 
interpretation of ^ in (3.38) as the "wave function" of the boundary state. 

The "crossing equations" (3.33), (3.35) and the expression (3.38) for the boundary 
state are the main results of this Section. We feel that the simple universal form (3.38) of 
the boundary state is not accidental. However at the moment we do not have a satisfactory 
understanding of its profound meaning and we can not offer anything better than the direct 
derivation through (3.31). 

4'. Now we turn to the "boundary bootstrap conditions". There are two sorts of these: 
the bootstrap conditions describing the boundary scattering of the "bound-state" particles 
and the conditions related to possible existence of the "boundary bound states" . 

If the particle A c can be interpreted as the bound state of A a A b (i.e. the pole at 9 = 
iu c ab , with u° ab satisfying (3.11)), the boundary S-matrix elements Rc(9) can be obtained 
by taking the appropriate residue in the bound-state pole of the two-particle boundary 
S-matrix R a \{9\, #2). This way one gets the equation 

ffR d M = fZ iai R%(0 + iu b ad )S b b f a2 (29 + iu b ad - tu a bd )R b b2 (9 - iu a bd ) (3.40) 

(u = in — u) which has the diagrammatic representation shown in Fig. 13. If the particles 
A a and A b have equal masses one can expect the appearence of a pole of R b i (9) at 9 = 

c 

T ~~ = ^ab ~ T due t° the diagram in Fig. 14. The corresponding residue can be 
written as 

• fab c 

K"\9)^ 1 -/^, (3.41) 
2 9- iu c ab 

where the "three-particle couplings" / are the same as in (3.12) but g° are new constants 
describing the "couplings" of the particles A c to the boundary (Fig. 14). More precisely, the 
nonzero value of g c indicates that the boundary state | B) contains a separate contribution 
of the zero-momentum particle A c , i.e. 

1 f°° 

I B) = M{1 + g c A c (0) + - / d9K ab (9)A a (-9)A b (9) + ...) | 0). (3.42) 

2 J -00 

Let us stress again that in our previous analysis of the boundary state we have ignored 
this possibility, i.e. strictly speaking validity of the equation (3.38) above is limited to the 
case when g° = for all the particles A c in the theory. It is easy to show that 
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[g c A c (0),K(9)} = 0, 



(3.43) 



where K{6) is the bilinear operator (3.39), and so in the general case one can look for 
the "boundary wave function" in the form ^[g c A c (0), K{6)\. The commutativity (3.43) 
follows from the relation 



g c 'K a ' b (9)S c c ? a ,(9) = g c 'K ab '(9)S$ c ,(9) (3.44) 

(Fig. 15) which is easily obtained if one considers the limit 9 — > ^ — \u c aihi in (3.25) and 
takes into account (3.13). It is also possible to show that if g Cl 7 g C2 ^ the amplitude 
K ClC2 (9) has the pole at 9 = with the residue 

K^(9)^- l -^p; (3.45) 

This pole term is illustrated by the diagram in Fig. 16. 

Let us illustrate these bootstrap conditions with an example of the so-called "Lie- 
Yang field theory" . This field theory describes the scaling limit of Ising Model with purely 
imaginary external field, the critical point being the "Lie- Yang edge singularity" (see [18]). 
It is a massive field theory which can be obtained by perturbing the c = — ^ minimal 
CFT by its only nontrivial primary field <p = $(i,2)? which has a conformal dimension 
A = — |. Mussardo and Cardyfll] have shown that the "bulk" theory is integrable; they 
have also found the corresponding factorizable "bulk" S-matrix. The theory contains only 
one species of particles, A, and the "bulk" S-matrix is described by 



A(9 1 )A(9 2 ) = S(9 1 - 9 2 )A(9 2 )A(9 1 ) (3.46) 

with 

, sinh 9 + i sin 

S(9) = 3.47 

V ' sinh6>-isin^ V ' 

The pole of S(9) at 9 = ( which has negative residue thus making non-unitarity of this 
theory manifest) corresponds to the same particle A appearing as the U AA bound state" . 
Here we do not attempt to analyze the possible integrable boundary conditions in this 
theory; we just assume that such ones do exist. Then the factorizable boundary S-matrix 
is described by (3.46) and 

A(9)B = R(6)A(-6)B, (3.48) 

where the boundary scattering amplitude R has to satisfy (3.26) and (3.35) (the boundary 
Yang-Baxter equation (3.25) is satisfied identically), i.e. 

R(9)R(-9) = 1- K{9) = S(29)K(-9); K{9) = R{— -9). (3.49) 

As the particle A appears as the "bound state" , the equation (3.40) has to be imposed as 
well, 
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R(0) = S(2B)R(6+ % -£)R(0-j). (3.50) 



There are two "minimal" solutions to (3.49), (3.50) 6 , 

sinh(f + f) sinh(f + sinh(| - f) 

R m (9) = ^| ^ ^| ^ i| 3.51 

(1)lJ sinh(|-f)sinh(|-^)sinh(f + f)' 1 ' 

R (9) - Sinh( ^ - ^ Sinh( ^ - ^ } sinh( * ' ^ f3 52) 

{2)y ) sinh(§ + f ) sinh(§ + ) sinh(§ + f ) l ' ) 

(of course R^ and _R( 2 ) are related through the "CDD factor" (3.27)). These two solutions 
have rather different physical properties. R^(9) exhibits a pole at 9 = ^ associated with 
the diagram in Fig. 14. This means that for this solution the boundary state | -B)(i) contains 
the single-particle contribution 

| fl) ( i) = (l + <7(i)A(0) + ...) | 0) (3.53) 
with non-zero amplitude g^iy Explicit calculation gives 



=2^2^3-3. (3.54) 
Accordingly, the amplitude R^(9) has a pole at 9 = with the residue 

R W m = (3.55) 

The solution R^ does not have any poles in the physical strip, i.e. g^ 2 ) = 0. 

In the general case the poles (3.41) and (3.45) do not exhaust all possible singularities 
the amplitudes R b a {9) may have in the "physical strip" < Im9 < tt. Even with given 
boundary action (2.8)((2.9)) the boundary can exist in several stable states | cx)b] ol = 
0, 1, ub — 1 (let us stress that all these states belong to Hb)- These states are eigenstates 
of the operators (2.16) 

H s | a) B = | a) B , (3.56) 

with some eigenvalues h^; in particular e a = ha^ is the energy of the state | ol)b- We 
assume that eo is the smallest of e a , i.e. as before the state | 0)s is the ground state of 
Hb • If e a — eo < minm a the state | a) b is stable just because its decay is forbidden 
energetically; higher IM H s could ensure stability of | a) b even if e a — eo > minm a . In 
any case, the states | a) b must show up as the virtual states in the boundary scattering 
processes. The amplitudes R b a (9) can exhibit the poles 



6 We call "minimal" the solution which has "minimal" set of zeroes and poles in the 
physical strip < Im9 < 7r, i.e. one can not reduce the total number of zeroes and poles 
in this strip without violating (3.49), (3.50). 

7 Here we assume that the ground state is not degenerate. 
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R b (9) ~ - / a ° ya , (3.57) 
where g^ are "boundary-particle couplings" and Vq u satisfy 

e + m a cos(i>o a ) = e Q , (3.58) 

see Fig. 17. Thus the states | q)b with e a > eo can be interpreted as the "boundary bound 
states" , e a — eo being the binding energy. 

In this situation it is natural to generalize (3.22) to 

A a (6)B a = R b J 3 a (8)A b (-6)B l3 , (3.59) 

where B a are formal "boundary creating operators" analogous to (3.21) and R b £ a {6) are 
amplitudes of the scattering process involving the change of the state of boundary, as 
shown in Fig. 18. Again, the amplitude R b £y{9) does not vanish only if m a = m,b and 
e a = ep, as it follows from (3.56). The equations (3.25), (3.26), (3.33), (3.35), (3.40) above 
can be generalized in an obvious way to include these amplitudes. Here we quote only the 
"boundary bound state bootstrap equation" 

g2a4?W) = diAr (* - ™L)Rt:Z(e)s:i' 2 (e + ivSa) (3.60) 

shown in Fig. 19; here the parameter v@ a satisfies the equation 

e a + m a cos(t>f Q ) = ep, (3.61) 

analogous to (3.38), and g^ a are the "particle-boundary coupling constants" (Fig. 20) en- 
tering the residues 

Z u lUaa 

Note that (3.44) and (3.45) can be considered as particular cases a = ,9 = 7 = 0of (3.60) 
and (3.62), resp., if we think of the ground state | 0)b as the "bound state" of some particle 
A c (with g c ^ 0) and itself. 

We have implicitely assumed in the above discussion that the ground state | 0)s is 
non-degenerate. This assumption is taken for the sake of simplicity only; it is not nesessary 
either from a physical point of view or for mathematical consistency. It is straightforward 
to incorporate the possibility of degenerate ground states into the above picture. In the 
next two Sections, where we consider examples of the factorizable boundary scattering 
theory, we will encounter this interesting possibility. 



4.ISING MODEL. 

As is known, the scaling limit of the Ising Model with zero external field is described 
by the free Majorana fermion field theory 
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A= J dydx ciff(iP,iP), (4.1) 

where 

a F F(ip, V>) = ipdzij> - i'dztf) + mtfjifj. (4.2) 

Here (z, z) are complex coordinates and m ~ T c — T. The field theories corresponding 
to the high-temperature (T — T c — > + ) and the low temperature (T — T c — > 0~) phases 
of the Ising Model are equivalent (they are related through the duality transformation 
ijj — > V> — > — V>, which changes the sign of m in (4.1)). Here we will assume that m > 
and interpret this field theory as the low-temperature phase. In this phase there are two 
degenerate ground states, | 0, ±), so that the corresponding expectation values of the spin 
field a(x) are {cr(x))± = ±cr, where a ~ ms is the spontaneous magnetization. 

The bulk theory (4.1) contains one sort of particles - the free fermion A with the mass 
m. Intuitively, this particle can be understood as a "kink" (or "domain wall") separating 
domains of opposite magnetization. The corresponding particle creation operator A' (9) 
(we denote it here A^ to comply with the conventional notations) can be defined through 
the decomposition 

/oo 
d6[u;e%A(6)e imxsinhe+imtcoshe + o e i A* (6)e- imxainhe - imtcoshe ]; 
-oo 

/oo 
d6[u)e-i A(e)e imxsinh0+imtcosh0 + ue~% A^ (e)e- %mxsinh0 - imt cosh0 ], (4.3) 
-oo 

where t = iy and u> = exp(^); uj = exp(— ^). The operators A(9), A< (9) satisfy canoni- 
cal anticommutation relations {A{9), A\9')} = 5(9 - 9'); {A{0),A{9')} = {A^ (6) , A\9')} 
= 0. The last of these relations, 

A\9)A\9') = -A\9')A\9), (4.4) 
has the same meaning as (3.3) so that the free-fermion two particle S-matrix is [19] 

S = -l. (4.5) 

Let us consider now this field theory in the half-plane x < 0, with the boundary at 
x = 0. Assuming that the boundary conditions are chosen to be integrable, we can define 
the boundary scattering amplitude R(9), 

A\9)B = R{9)A\-9)B. (4.6) 
As is discussed in Sect. 3, this amplitude has to satisfy (3.26), (3.33), (3.35), i.e. 

R(9)R(-9) = 1; K(9) = -K(-6); 
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K(e) = R( l ^-e). (4.7) 

We consider first two simplest boundary conditions - the "free" and the "fixed" ones 8 . 

a). "Fixed" boundary condition. In the microscopic theory this boundary condition 
corresponds to fixing the boundary spins to be, say, +1. Obviously, this boundary condition 
removes the ground state degeneracy. In terms of the fermion fields ip, ip the "fixed" 
boundary condition can be written as 

(^ + V0.=o = O. (4.8) 

In presence of the boundary, the fields tfj^tf) still enjoy the decomposition (4.3), although 
the operators A, A' (now acting in He) are not all independant but satisfy the relations 

(ue% + ue~%)A(9) = -(u)e% + uje~^)A{-9); 

(uei +uje~^)A^(9) = -(we' + ue~^)A^(-9) (4.9) 
which follow from (4.8). From the last of (4.9) we find 

Z7T 

R fiX ed(B) = itanh(—--). (4.10) 

Alternatively, one could use another hamiltonian picture, with r = — ix interpreted as the 
time. In this picture the fields x = ^V*? X = ^ en j°y the same decomposition (4.3) with the 
substitution x — > y,t — > r and operators A(9), A^(9) acting in the Hilbert space Tt of the 
bulk theory. The boundary condition (4.6) appears as the initial condition (x + zx) T= o = 
which is understood as the equation 

( X + ix)r=o I B fixed ) = (4.11) 
for the corresponding boundary state | Bf ixed ). In terms of A, A^ this equation reads 

[cosh(9/ 2) A(8) +zsinh(0/2)At(-0)] | B flxed ) = 0, (4.12) 
and hence the boundary state can be written as 

I B flxed ) = exp{- / d9K fixed {9)A\-e)A\e)} \ 0) (4.13) 

with 

Q 

Kf lxed {9) = ztanh- (4.14) 

and | 0) =| 0, +). Although (4.13) with | 0) =| 0, — ) solves (4.12) as well, it is easy to see 
that in the infinite system this state does not contribute to | Bf ixed ) (in a large system, 



8 In the conformal field theory of the Ising model (which describes the case m = 0), 
these are just the two possible conformal boundary conditions; they are analyzed in [7,8]. 

20 



finite in the y-direction, —L/2 < y < L/2, its contribution is suppressed as exp(— mL)). 
Note that (4.10), (4.14) satisfy (4.7). 

b) . "Free" boundary condition. In the microscopic theory one imposes no restrictions 
on the boundary spins. Correspondingly, the ground state is still two-fold degenerate. 
Again, in terms of the fermions ip, ifj this boundary condition has a very simple form 

(V>-VO*=o = 0. (4.15) 
The same computation as in the previous case gives 

Rfree(0) = -icoth(---). (4.16) 

Note that the corresponding boundary state amplitude 

K free (0) = -icoth d - (4.17) 

exhibits a pole at 6 = (this pole is related to the existence of the zero-energy mode 
ip = -0 ~ expmi which satisfies (4.15)), which indicates that the boundary state | Bf ree ) 
contains the contribution of a zero-momentum one-particle state, 

| B free ) = (1 + ^(0) + ...) | 0). (4.18) 

Of course, this feature is easily understood in physical terms. Semi-infinite Ising Model at 
T, T c with free boundary condition admits a particular equilibrum state, characterized by 
the asymptotic conditions (a(x,y)) — > +a as y — > +oo and (a(x,y)) — > —a as y — > — oo. 
This state contains an infinitely long (fluctuating) "domain wall" attached to the boundary, 
which separates two domains of opposite magnetization, as shown in Fig. 21. This "domain 
wall" configuration is interpreted as a zero-momentum particle emitted by the boundary 
state. It is not difficult to check that the boundary state 

| B free ) = (1 + At(0))exp{- / K free {9)A\-9)A\e)} | 0) (4.19) 
satisfies the boundary state equation 

J d6f(6)[smh(6/2)A(6)-icosh(6/2)A\-6)] | B free ) = (4.20) 

(f(0) is an arbitrary smooth function) which follows from the "initial condition" (x — 

ix)r=o = 0. 

The simple boundary conditions above can be obtained as the two limiting cases of 
the more general integrable boundary condition which we describe below. 

c) . "Boundary magnetic field". This boundary condition is obtained by introducing a 
nonzero external field h ( "boundary magmetic field" ) which couples only to the boundary 
spins. Obviously, h = correspond to the "free" case above while in the limit h — > oo one 
recovers the "fixed" boundary condition. The boundary magnetic field can be considered 
as the perturbation of the "free" boundary condition, 
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/oo 
dya B {y), (4.21) 
-oo 

where (Jb{u) is the "boundary spin operator" [7,8] 9 . The action Af ree is 

/oo />0 foo 

dy dxa FF {^^) + 1/2 / %[(#) I=0 + H- (4.22) 
-oo J — oo J — CO 

Here d = ^a, and a(y) is an additional (fermionic) boundary degree of freedom which is 
introduced to describe the ground-state degeneracy. The quantum operator a associated 
with the boundary field a(y) anticommutes with ip, ?/> and satisfies 

a 2 = 1, (4.23) 
so that for the free boundary condition one has 

a\0,±) B =\0,T) B . (4.24) 

The equation (4.15) follows directly from the action (4.22). The boundary spin operator 
cb(i/) is analyzed (along with other boundary operators) in [7,8]. It is identified with the 
degenerate primary boundary field ^(1,3) °f dimension A = 1/2. In the Lagrangian picture 
described above it can be written as 

a B (y) = \{i> + $)x=o(y)a(y), (4.25) 

so that the dimension of h in (4.21) is [mass] 1 / 2 . 

One obtains, from (4.21) and (4.25), the boundary condition for the fermi fields ip, 

d h 2 
^(V>-V0*=o= y(V> + V0*=o- (4-26) 

With this, it is straightforward to compute the boundary scattering amplitude for (4.21), 

R h (6) = ztanh(z7r/4 - 9/2) - ~ \ ^ (4.27) 

k + i smh 6 

where 

h 2 

k = 1 . 4.28 

2m v ' 

Note that for h > the amplitude Rh(0) does not have a pole at 6 = in/ 2 and hence 
there is no one-particle (and any odd-number particle) contributions to the corresponding 
boundary state | Bh); this state has therefore the general form (3.38) with | 0) =| 0, ±), 
for h = ±|/i|, resp., and with 



9 This perturbation generates a "flow" from the free boundary condition down to the 
fixed one. In the case of conformal "bulk" theory (m = 0) this flow was studied in [20] . 
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K h {9) =ztanh(^/2) /t + C ° S ^ . (4.29) 
k — cosh 9 

Again, this is not very surprising. The nonzero boundary magnetic field removes the 
ground-state degeneracy; for h sufficiently small, the two-fold degenerate ground state 
| 0, ±}b of the "free" boundary splits into two non-degenerate states | 0)b and | 1)b, 
where | 1)b can be interpreted as the boundary bound state. For < h 2 < 2m we can 
parametrize (4.28) as 

k = cos v. (4.30) 

with real v, < v < n/2. In this domain the amplitude (4.27) exhibits a pole in the physical 
strip at 9 = i(iv/2 — v) which is associated with the state | 1)^. Physical interpretation 
of this boundary bound state is simple. For h ^ the equilibrium "ground state" which 
minimizes the free energy features an asymptotic behaviour (a(x,y))o — > +sign(h)a as 
x — > — oo. However, for h sufficiently small, there exists another stable equilibrum state 
with (a(x, y))i — > —sign(h)a as x — > — oo. Although its free energy is higher by a positive 
boundary term, it is indeed stable as in the infinite system there is no finite kinetics of decay 
(in a system finite in the y direction the decay probability is suppressed as exp(— mL cos v)). 
Clearly, the energies eo and e\ of the states | 0)b and | 1)b have the meaning of specific 
(per unit boundary length) boundary free energies of these two equilibrum states. From 
(4.27) we have 

e\ — eo = msinv. (4-31) 

It is also clear that the expectation values can be obtained by analytic continuation 
of (...)o from h to — h. So, in this very domain \h\ < \[2~m "boundary histeresis" [21] can 
be observed. It is possible to show that these expectation values can be computed as the 
matrix elements 



<■■■>- -^mr-'-w^' (432) 

where the last expression contains the "excited boundary state" 

I B' h ) = exp(^ I d9K h {9)A\-9)A\e)) \ 0') (4.33) 

with the integration taken along contour C shown in Fig. 2, and | 0') =| 0, =f) for h = ±|/i|. 

As h approaches the "critical" value h c = \j2m (i.e. v approaches 7r/2), the "boundary 
bound state" becomes weakly bound; e\ — eo — m = — 2msin 2 (7r/4 — v/2) — > 0, and its 
effective size £ = (ei — eo — m) _1 diverges as (h c — h)~ 2 . Correspondingly, the equilibrium 
state developes large boundary fluctuations which propagate deep into the bulk, 

(cr(x, y))i + a ~ exp(x/^) as x — > — oo. (4.34) 
At h = h c the boundary condition (4.26) reduces to 
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(4.35) 



and the boundary scattering amplitude can be written as 

R h =h c (9) = -i tanh(z7r/4 - 9/2) (4.36) 

which differs only by a sign from (4.10). This difference is very significant though; (4.36) 
admits the state A^(0) \0)b - remnant of | 1)b- 

For 2m < h 2 < 4m we can still parametrize k as in (4.30), with real v. However, 
now 7r/2 < v < 7r. The pole of Rh at 9 = i(n/2 — v) leaves the physical strip; there is no 
"boundary bound state" in this domain, and the equilibrum state is unique. At h = 2y/m 
the expression (4.27) simplifies as 

R h =2^{0) = -ztanh 3 (z7r/4 - 9/2), , (4.37) 

i.e. the boundary-state amplitude K h=2y /tE possesses a third-order zero at 9 = 0. This 
seems to affect the x — > — oo asymptotic of (cr(x, y)), but at the moment we do not have a 
clear physical interpretation of this phenomenon. 

Finally, at h 2 > 4m one can write k = — cosh 9q with 9q real and positive. In this 
regime the amplitude Rh shows two "resonance" poles at 9 = — ire/ 2 ± 9q. As h grows 
these poles depart to infinity and ultimately Roo reduces to (4.10). 

Clearly, the behavior described above agrees with the expected "flow" from the "free" 
boundary condition down to the "fixed" one[20]. This example is particularly simple 
because the described boundary field theory is not only integrable but free - both the 
bulk and the boundary equations of motion for ip, ij) are linear. It is easy to show that all 
solutions of (4.7) correspond to free-field boundary conditions, generalizing (4.26) in an 
obvious way. 



5. BOUNDARY SINE-GORDON MODEL 

In this Section we consider the integrable boundary theory for the sine-Gordon model. 
The bulk sine-Gordon field theory is described by the euclidean action (2.1) with 

1 m 2 

a(<p, d M v?) = g ( d ^) 2 - cos(/3(p) (5.1) 

where ip(x,y) is a scalar field and (3 is a dimensionless coupling constant. This model is 
well known to be integrable both as classical and quantum field theories (see e.g. [22]. In 
the quantum theory, the discrete symmetry tp-^Lp+^-N,NE Zis spontaneously broken 

at f3 2 < 87r[23]; in this domain the theory is massive and its particle spectrum consists of 
a soliton-antisoliton pair (A, A) (with equal masses) and a number (which depends on f3) 
of neutral particles ("quantum breathers") B n: n = 1,2, ... < A, where 
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The soliton (antisoliton) carries a positive (negative) unit of "topological charge" 

P f + °° d (3 
q = — dx—(p(x,y) = —[cp(+oo,y)-cp(-oo,y)]. (5.3) 

27T J-oo OX ZTT 

The charge conjugation C : A <-> A is related to <p <-> — ip symmetry of (5.1). The particles 
B n are neutral (they are interpreted as the soliton-antisoliton bound states), B n with even 
(odd) n being C -even (C -odd). The masses of B n are 

m n = 2M s sin( — ); n = l,2,... < A, (5.4) 

where M s is the soliton mass. 

Factorizable scattering of solitons is described by the commutation relations 

A{6)A{6') = a(0 - 0')A(e')A(0), 1(0)1(0') = a(9 - 9')A(9')A(9), 

A(0)A(0') = b(0 - 0')A(0')A(0) + c(0 - 6')A(e')A(9), (5.5) 

where A(0) and A(0) are soliton and antisoliton creation operators and the two-particle 
scattering amplitudes a, 6, c are given by 

a(9) = sin(A(7r — u))p(u), 

b(9) = sm(Xu)p(u), (5.6) 
c(9) = sin(A7r)p( , u), 

where u = —i9 and 

P (u) = -Ir(A)r(i - XuMT(i - a + xu/,) f[ 

pr v r(2/A - \u/n)r(l + 2/A - Xu/n) 

1{U) r((2/ + 1)A - Au/7r)r(l + (2/ - 1)A - Au/tt) " 1 J 

The amplitudes of A£? n and B n B m scatterings can be obtained from these with the use of 
(3.13), they can be found in [4]. 

Let us consider now this field theory in presence of the boundary at x = 0, described 
by the action (2.8). The first question which arises here is how to choose the "boundary 
action" 6(</?s, d/dy^s) in (2.8) in order to preserve integrability. General classification of 
integrable boundary actions seems to be an interesting open problem. Here we just claim 
that the boundary sine-Gordon theory (2.8), (5.1) with 

b(<p) = -Mcos(^(<p-<p )), (5.8) 

where M and po are free parameters, is integrable. This claim is supported by an explicit 
computation of the first nontrivial integral of motion in the Appendix, where we restrict 
attention to the classical case; we believe this computation can be extended to the quantum 
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theory along the lines explained in Sect. 2. We want to describe the factorizable boundary 
scattering theory associated with (5.8). 

The soliton (antisoliton) boundary scattering amplitudes can be encoded in the com- 
mutation relations (see Section 2) 

A(9)B = P+(6)A(6)B + Q+(6)A(6)B; 

A(6)B = P_(6)A(6)B + Q_(6)A(6)B. (5.9) 

Here P+, Q+ (P_, Q-) are the amplitudes of soliton (antisoliton) one-particle boundary 
scattering processes shown in Fig. 23. Some comments are in order. Exept for the case 
M = oo (see below), the boundary value <p(x = 0,y) is not fixed in the boundary field 
theory (5.8)and hence the topological charge 

(3 f° d 

q= 2^]_J X d- X ^ v) (5 - 10) 

is not conserved. Therefore the boundary scatterings are allowed to violate this charge 
conservation by an even number. In particular, an incoming soliton can go away as an 
antisoliton and vice versa. The amplitudes Q± are introduced to describe these processes. 
The excceptional case is M = oo ( "fixed" boundary condition) ; in this case we must have 
Q± = 0. Also, at <po 7^ 0(mod27r / (3) (and M ^ 0) the boundary interaction violates 
charge-conjugation symmetry; that is why in the general case the amplitudes P+, Q+ and 
P_, Q- are expected to be different. 

With (5.9), the boundary Yang-Baxter equation (3.25) leads to six distinct functional 
equations, 



Q+(6')c(6 + e')Q-{9)a(e - 6') = Q_(0')c(0 + 6')Q + (e)a(e - 9'), (5.11a) 

P-(6')b(6 + 6')P+(6)c(6 - 9') + P+(9')c(9 + 9')P+(9)b(9 - 9') = 
P-(9')c(9 + 9')P_(9)b(9 - 9') + P+(9')b(9 + 0')P_(0)c(0 - 0'), (5.116) 

P + (0Oa(0+0OQ + (0)c(0-0O+Q+(0OK^O^+(^)K^-^)+Q+(^)c(0+0 / )^-Wc(0-0 / ) = 
Q+(0>(0 + 0')P+(0)a(0 - 9') + P_(0')c(0 + 9')Q + (9)a(9 - 9') (5.11c) 

P + (0Oa(0+0OQ + (0)6(0-0O+Q+(0OK^O^+(^c(0-0O+Q+(^Oc(0+0O^-W6(0-0 / ) = 

P+(9')b(9 + 9')Q+(9)a(9 - 0'), (5.11d) 

(the other two are obtained from (5.11c), (5. lid) by + <-> — substitution). The solution 
to these equations is found in a recent paper [24] 2 ; it reads 

P + (0) = cos(£ + Xu)R(u); 
P_(0) = cos(£-A-u)P(w); 

2 We learned about [24] after we had found this solution independently. 
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Q+(6) = ^-sm(2\u)R(u); 

Q_(6) = ^-sm(2\u)R(u), (5.12) 

where again u = —i9; £, k± are free parameters and R(u) is an arbitrary function 3 . 

The free parameters in (5.12) have to be related to the parameters of the boundary 
action. The solution (5.12) seems to contain one parameter more than the action (5.8). It 
is easy to see however that one of the two parameters k± in (5.12) can be removed by an 
appropriate gauge transformation 

A{6) -> e lQ A(6>), A{6) -> e - za A(6) (5.13) 

which leaves all the charge-conserving amplitudes unchanged and transforms the ampli- 
tudes Q± as 

Q+(0) -> e- 2ta Q + {6)- Q_(0) -> e +2ia Q-{6). (5.14) 

Considered from the Lagrangian point of view, this transformation amounts to adding a 
total derivative term to the boundary action density (5.8), 

%)-%) + ^^- (5-15) 

We use this freedom to set k+ = k- = k. This leaves us with two parameters, £ and k, 
which correspond to the two-parameters M and <po in (5.8). 

The function R(u) in (5.12) can be determined with the use of the "boundary unitarity" 
and the "boundary cross-unitarity" equations (3.26) and (3.35). These equations reduce 
to 

R(u)R(-u) = [cos 2 (0 - sin 2 (Aw) - k 2 sin 2 (Aw) cos 2 (Au)] _1 ; (5.16) 

K(u) = sm\(n + 2u)p(2u)K(-u); K(u) = R{tt/2 - u). (5.17) 

It is not difficult to solve these equations. In fact, it is convenient to solve (5.16) and (5.17) 
"separately" . Let us factorize the function R(u) as 

R(u) = R (u)R!(u); K(u) = Ko^K^u), (5.18) 
where K (u) = R (n/2 — u), Ki(u) = R\{tt/2 — u) so that these factors satisfy 

Ro(u)Rq(-u) = 1; K (u) = sinA(7r + 2u)p(2u)K (-u) (5.19) 

and 



i*i(u)fli(-u) = [cos 2 (0 - (1 + k 2 ) sin 2 (Aw) + k 2 sin 4 (Au)] _1 ; Kx{u) = K^-u). (5.20) 
3 This is the general solution for generic A. For integer A there are additional solutions. 
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The equations (5.19) do not contain the boundary parameters £ and k; the solution can 
be written as 



R (u) = F (u)/F (-u); 

w . r(l - 2Xu/n) 
Fo{u) = T(X-2XuM X 

yr T(4Xk - 2Xu/n)T(l + 4Xk - 2Xu/>tr)T(X(4k + 1))T(1 + X(4k - 1)) 

r(A(4fc + 1) - 2Aw/7r)r(l + X(4k - 1) - 2Aw/7r)r(l + 4Xk)T(4Xk) ( ' ^ 

The factor R contains the poles in the "physical strip" < u < tv/2 located at u = u n = 
mv/2X;n = 1,2, ... < A. Appearence of these poles is very well expected. In the generic 
case, the boundary state | B) associated with the boundary condition (5.8) is expected 
to contain the contributions of the zero- momentum particles B n , which leads to the poles 
(3.41). So, these poles of R correspond to the diagrams in Fig. 24. Note that these poles 
should not appear in the amplitudes Q± and the factor sin(2Aw) in (5.12) takes care of 
this. 

The equation (5.20) contains all the information about the boundary condition (i.e. 
the parameters £ and k). Its solution can be written as 

R x (u) = -^—a(j], u)a{id, u) (5.22) 
cos 4 

where 



_ U(x, tt/2 - u)U(-x, 7r/2 - u)U(x, -tt/2 + u)IL(-x, -tt/2 + u) 
a[X ' U) ~ n 2 (x,7r/2)n 2 (-a:,7r/2) 

_ TT r (V 2 + (gj + V 2 ) A + x / n ~ ^u/n)T(l/2 + (21 + 3/2)A + x/n) 
~ to 17(1/2 + (2/ + 3/2)A + x/7r ~ Aw A)r(V 2 + (21 + 1/2)A + x/tt) 

solves 



a(x, u)a(x, —u) = [cos(x + Xu) cos(x — Xu)] l ; a(x, n/2 — u) = a(x, tt/2 + it), (5.24) 
and the parameters r\ and $ are determined through the equations 

cos (77) cosh (#) = — ^ cos £; cos 2 (?7) + cosh 2 (z?) = 1 + -j^. (5.25) 

The above boundary S-matrix has a very rich structure. The factor a(i"d, u) in (5.22) 
brings in an infinite set of singularities at complex values of u; some of these complex 
poles probably can be interpreted as resonance states of the boundary. The factor 17(77, u ) 
exhibits infinitely many poles at real u; depending on the parameter 77, some of these poles 
can occur in the "physical strip" < u < n/2 thus giving rise to the boundary bound 
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states. We did not carry out a complete analysis of the possible boundary phenomena 
described by this S-matrix. In order to do this one must find the an exact relation between 
the parameters M and (po in the boundary action and £ and k (or rj and $) in the S-matrix. 
A brief analysis shows that in the general case this relation is very complicated. We leave 
this problem for future studies. Here we analyze only two particular cases. 

a) . "Fixed" boundary condition <p(x,y) x =o = po can be obtained from (5.8) in the 
limit M = oo. In this case the topological charge is conserved and the amplitudes Q± in 
(5.9) must vanish. Hence this case corresponds to k = 0. We have two amplitudes, P+ and 
P_ which describe soliton scattering processes schematically shown in Fig. 25. For k = 0, 
the "resonance" factor a(ii}, u) disappears and the equation (5.22) simplifies as 

Ri(u) = ^—a(^u), (5.26) 
cos£ 

where £ is related in some way to ipo. Obviously, ipo = corresponds to £ = as at (po = 
the boundary theory respects C symmetry and we must have P+ = P_. At £ = the 
function Ri(u) does not exhibit any poles in the physical strip (and hence there are no 
boundary bound states), the nearest singularity being the double pole at u = —iv/2X. At 
£ 7^ (i.e. <po 7^ 0) this double pole splits into two simple poles at u = u± = — 7r/2A±£/A. 
Note that due to the factors cos(£ ± Xu) in (5.12) the pole 9 = iu + (9 = iu-) appears only 
in P + (P-). At £ > tv/2 the pole at 9 = iu + enters the physical strip. Appearence of such a 
"boundary bound state" is well expected. For < (p < n//3 the ground state | 0)s of the 
boundary sine-Gordon theory is characterized by the asymptotic behaviour (<p(x,y)) b — > 
as x — > — oo. Classically, there is another stable state with {(p(x,y)) — > 2tv/(3 asi^ — oo. 
If (po is small (compared to the quantum parameter /?), quantum fluctuations can destroy 
its stability. However, if </? is not too small, this state is expected to be stable in the 
quantum theory, too. At <po = ir/fl this state has the same energy as the ground state 
which hence becomes degenerate. In this case the amplitude P+(0) has to have a pole at 
9 = in/2 corresponding to the "emission" of a zero-momentum soliton by the associated 
boundary state, as is explained in Sect. 3. We can conclude that tpo = ir/f3 corresponds to 
the case u + = tv/2. Assuming a linear relation between £ and </3 , we find 

£ = J ¥>o. (5-27) 

Let us stress that the assumed linear relation (and hence (5.27)) is no more than a con- 
jecture. Also, the relation (5.27) is conjectured only for the "fixed" case M = oo; it is not 
expected to hold in the general case of finite M. 

b) . "Free" boundary condition corresponds to M = in (5.8). In this case the full 
action enjoys the symmetry cp — * — <p (C symmetry) and hence we must have P + = P- = 
Pfree and Q+ = Q- = Qfree, i-G. we have to choose £ = 0. Note that under this choice the 
factors cos(A#) in (5.12) cancel the poles at 9 = iu2 n +i of the factor (5.21). This agrees 
with the expectation that in the C -symmetric case the boundary state can not emit the 
C -odd particles B 2n +i S - In addition, in the "free" case one expects all the amplitudes 

3 In fact, it is easy to argue that in general the choice £ = corresponds to the C 
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(5.12) to have a pole at 6 = nr/2, for exactly the same reason which we discussed in the 
previous Section. To ensure this property one has to choose rj = §(A + 1) (and d = 0), i.e. 
k = [sin(7rA/2)] _1 . So in the "free" case the amplitudes (5.12) can be written as 



Pfree(0) = SUl(7r X/2)Rf ree (u) ; 



QfreeiO) = Sm(Xu)Rf ree (u) 



(5.28) 



where 




cos(Aw) 



Ro(u)<r(4ir 2 //3 2 ,u)a(0, u). 



(5.29) 



sin(A7r/2) 



6. DISCUSSION 

Obviously, there is much room for further development. 

More detailed analysis of "integrable boundary conditions" is needed. In Section 
2 we assumed that the boundary action b in (2.8) depends only on boundary values 
<Pb(h), d/dy(pB(y) of the "bulk" field <p(x, y). In general, the boundary action can contain 
also specific boundary degrees of freedom. Clearly, this provides more possibilities for the 
integrable boundary field theories. 

Further study of possible solutions to Boundary Yang-Baxter equations ((3.25) is of 
much interest. A new class of solutions is found in a recent paper [24]. 

Many integrable "bulk" theories exhibit nonlocal integrals of motion of fractional spin 
[25-27], along with the local ones. Analysis of these nonlocal integrals provides valuable 
information about the S-matrix. We expect that the nonlocal integrals of motion exist in 
some integrable boundary field theories. A natural place to look for the nonlocal IM is 
boundary sine-Gordon model considered in Section 5. As in the case of bulk sine-Gordon 
theory, where nonlocal IM generate the so-called "affine quantum group symmetry" [27], 
appropriate modifications of these IM could help to find an exact relation between the 
parameters of the action (M and <po in (5.8)) and of the S-matrix (£, k in (5.12))(which 
we found very difficult to guess). 

The most interesting problem is how to extract any off-shell data from the S-matrix. 
In the bulk theory, two approaches have proven to be the most successfull. One is the 
"formfactor bootstrap" proposed by Smirnov[27]. It would be interesting to generalize the 
formfactor bootstrap equations to the case of integrable boundary field theory, both for 
bulk operators in presence of boundary and for boundary operators. Another approach is 
known as "thermodinamic Bethe anzats" (TBA)[29]. At the present stage of development, 
it is capable of providing the ground state energy of a finite-size system (with spatial 
coordinate compactified on a circle), once the S-matrix is known. This approach can 
be generalized in different ways to incorporate the boundary effects. In particular, the 
expression (3.38) for the boundary state allows one to find the analog of TBA equations 
for spatially finite systems with integrable boundaries. We intend to study these equations 
elsewhere. 

-symmetric boundary action (5.8) with (po = 0. Even in this case the relation between M 
and k seems to be rather complicated. 
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It was understood recently that many interesting "massless flows" in the bulk the- 
ory can be described in terms of "massless factorizable S-matrices", through TBA equa- 
tions [30]. This approach can be modified to incorporate "massless boundary S-matrices" 
[31], thus providing the possibility of a similar description of the "boundary flows" [32]. 

After this work was finished we received a recent paper [31] where some results of 
Section 3 (notably, the "boundary bootstrap equation" (3.40)) are obtained. However, the 
"boundary cross-unitarity equation" (3.35) (which we think to be the most significant of 
our results) is not considered there; that is why we decided that our paper is still worth 
publishing. 



APPENDIX 

Here we carry out a preliminary study of integrable boundary conditions in the sine- 
Gordon model (5.1). We restrict attention to the classical case (3 — > 0. In this limit it is 
convenient to work with the field (f>(x, y) = (3<fi(x, y) so that the "bulk" equation of motion 
takes the form 

d z d z (p = sm(p, (Al) 

where we set m 2 = 4. As is known [22], there are infinitely many conserved currents 
(T s+1 , 6 s _i) and (T s+1 , O s _i) with s = 1, 3, 5, which satisfy (2.5) in this "bulk" theory; 
the first two have the form 

T 2 = (a^) 2 ; e = -2cos(/>, (A2) 

T 4 = (d 2 z <j>) 2 - ^(d 2 </>) 4 ; 2 = («9^) 2 cos</> (A3) 

(T and are obtained from (A. 2) by substitution z — > z). The current (A. 2) is just the 
energy-momentum tensor, and (A. 3) leads to the first nontrivial IM. 

As is explained in Section 2, the current (A. 3) would lead to a nontrivial IM in presence 
of the boundary, too, provided the boundary condition at x = is chosen in such a way 
that 

[T 4 + e 2 - T 4 - e 2 ] x=0 = -^-03, (A4) 

dy 

where 9s is some local boundary field. Let us consider the simplest case of the boundary 
action b independant of the derivative i.e. the total action of the form 

A= Jv{ dy dx[-(d^(x,y)) 2 -4cos(j)(x,y)]+ / dyV(<f>(x = 0, y))}. (A5) 

P J — oo J — oo " J — oo 

The boundary condition which follows from (A. 5) is 

[d x( p + V'((P)] x = = 0. (A6) 
31 



With this, the l.h.s of (A. 4) can be written (up to an overall numerical factor) as 

V'{j)){d y $f + 8V"(<f>)(d 2 yy <f>)(d y cf>) -[(V(cf>)f + 4y // (0)sin0 + 2y / (0)cos0](^0). (A.7) 
This reduces to a total d/dy derivative if the function V(4>) satisfies 

4V"(<P) + V(<P) = 0, (A8) 

that is 

K(0) = -Acos(^^), (A9) 
where A and </>o are constants. Returning to the original normalization, one gets (5.8). 
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